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ABSTRACT 

To  highlight  in  a  qualitative  manner  the  kinetics  of  an 
excimer  laser,  a  simple  computer  model  for  calculating  the 
species  populations  in  a  KrF  laser  cavity  is  developed;  sub- 
sequently a  computer  program  originally  developed  at  the 
Naval  Research  Laboratory  (NRL)  is  modified  to  calculate  the 
population  of  the  different  electronic  configurations  of 
excited  neon  present  in  a  XeF  laser.   'v^Hien  modified,  the  NRL 
program  accounts  for  135  reactions  and  requires  9.5  minutes 
of  CPU  time  using  the  IBM  360/67.   The  populations  obtained 
are  applied  to  the  calculation  of  the  index  of  refraction 
in  the  laser  cavity.   The  phase  shift  is  determined  per  unit 
length  for  the  two  laser  wavelengths;  one  laser  wavelength  is 
non-resonance  with  neon  and  the  other  is  at  resonance.   Neon 
is  the  dominant  specie  relative  to  concentration  and  within 
the  population  distribution  the  neon  ground  state  dominates 
by  a  factor  of  a  million.   The  calculations  show  that  the 
ground  state  neon  determines  the  index  of  refraction;  an 
exception  occurs  if  the  laser  wavelength  is  near  resonance 
to  any  of  the  transitions  in  the  3s  ■>  4p  array.   As  long  as 
the  laser  is  operated  away  from  the  3s  ■+  4p  resonant  wave- 
lengths, the  phase  shift  will  be  negligibly  small  resulting 
in  satisfactory  beam  quality;  the  preceding  statement  is 
valid  only  for  the  influence  of  neon. 
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I.   INTRODUCTION 

The  rare-gas  monohalides ,  such  as  XeF ,  are  simple  diatomic 
molecules  whose  properties  and  emission  spectra  were  essen- 
tially unknown  as  recently  as  five  years  ago.   Now  they  are 
the  active  media  for  gas  lasers  that  could  provide  overall 
electrical  efficiency  as  high  as  10  per  cent  in  various  wave- 
lengths,  as  depicted  in  Table  I. 

Table  I.   Rare-gas  Halide  Laser  Emission  Wavelengths 
Peak  wavelengths  of  most  intense  band  (nm) 

F 

CI 
Br 
I 

Has  not  demonstrated  laser  oscillation. 

Predissociates ,  hence  emission  is  weak  or  unobservable . 

The  overall  efficiency  as  explained  by  J .  J.  Ewing  and  C.  A 

Brau  [1]  is 

...a  product  of  a  quantum  efficiency,  an  upper 
laser  level  production  efficiency,  an  extraction 
efficiency,  an  efficiency  for  producing  the  initial 
excited  states  by  the  pump,  and  an  energy  coupling 
efficiency  which  describes  the  efficiency  with  which 
energy  gets  from  a  wall  socket  into  the  gas  medium. 
The  extraction  efficiency  is  comprised  of  both 
spectroscopic  and  kinetic  extraction  efficiencies 


Ne 

Ar 

Kr 

Xe 

107a 

193 

248 

351,353 

b 

170 

222 

308 

b 

166a 

,b 

206a 

282 

b 

b 

185a: 

,b 

252a 

depending  on  the  ratio  of  net  gain  to  loss  and  the 
ratio  of  rates  of  stimulated  emission  of  upper  laser 
levels  to  the  quenching  and  spontaneous  loss  of 
upper  levels.   Certain  rare  gas  halide  lasers  have 
potential  for  high  power  and  high  efficiency  because 
they  have  high  efficiency  in  each  of  the  above 
mentioned  elements.   Depending  on  the  application 
requirements,  one  may  accept  a  lower  efficiency  of 
one  of  these  elements  if  the  compromise  allows  for 
some  simplification  of  the  laser.   Typically  this 
trade  off  usually  involves  coupling  and  excited  state 
production  efficiencies . 


Rare-gas  monohalides  belong  to  a  larger  class  of  mole- 
cules broadly  known  as  excimers .   An  excimer  is  an  atomic 
or  molecular  aggregate,  composed  of  two  atoms  or  molecules, 
that  is  bound  in  its  excited  state  but  is  unstable  or 
slightly  bound  in  the  ground  state. 

The  term  "excimer"  is  standard  chemical  nomenclature 
for  dimers  which  are  bound  in  the  excited  state  and  free  in 
the  lower  state.   In  standard  nomenclature  dimer  means  two 
of  the  same  thing.   Thus,  Xe^*,  Ar«* ,  etc.  are  excimers. 
Species  that  are  heteronuclear  such  as  LiXe*  or  KrF*  also 
can  have  bound  excited  states  and  dissociative  lower  states 
The  accepted  chemical  nomenclature  for  such  species  is 
"exciplex"  short  for  excited  complex.   Unfortunately,  the 
laser  community  has  not  uniformly  utilized  good  scientific 
English  and  typically  calls  all  species  with  dissociative 
lower  states,  and  even  some  like  XeF  with  slightly  bound 
lower  states,  "excimers". 

The  excited  state  of  the  excimer  species  can  radiate 
in  a  broad  band  that  is  red  shifted  from  the  wavelength  of 
the  parent  atomic  excitation.   The  size  of  the  bandshift 
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depends  on  the  depth  of  the  excited  state  well  (typically 
l-5eV)  and  on  the  repulsion  of  the  lower  state  at  the  energy 
minimum  of  the  excimer.   In  Figure  1  are  schematic  potential 
energy  curves  for  an  excimer. 


o 
w 


AB*   Moleculation 


ther  excited  states 
A+B* 


Excimer  upper  level 


Atomic 
excitation 


Excimer  emission 


A+B 


Continuum 
bandwidth 


Weak  van  der  Walls  binding 


INTERNUCLEAR  SEPARATION 


Figure  1.   Potential  Energy  of  Excimers  (Reproduced  from 
Ewing  [ 2 ] . 


Excimer  emission  spectra  are  useful  for  lasers  because 

the  lower  level  rapidly  dissociates.   Since  the  dissociation 

-12 
time,  about  10     second,  is  much  less  than  upper  level 

-9      -6 
radiative  lifetimes,  10    to  10    second,  population  inver- 
sions can  be  readily  produced  and  maintained. 

The  excimer  laser  or  exciplex  laser  offers  features 
which  are  essential  for  military  applications  [2-4].   These 
features  are  as  follows:   first,  the  potential  for  high 
efficiency;  second,  scalability  to  large  volume  and  high 
power;  and,  third,  transmission  through  the  atmosphere 
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using  either  directly-generated  photons  or  frequency-con- 
verted photons  [5].   As  a  result  of  the  attractive  character- 
istics, considerable  effort  has  been  and  is  being  directed 
toward  identification  of  key  technical  issues  and  establish- 
ment of  feasible  solutions  to  the  problems  identified  [6-11]. 
In  addition  to  research  directed  toward  the  laser  device 
itself,  work  is  also  being  funded  on  related  and  relevant 
aspects,  including  e-beam  technology,  UV  optics,  and  time 
dependent  laser  aerodynamics.   For  example,  see  Bennett, 
et  al .  [12]  for  work  on  UV  optics. 

Systems  studies  have  been  conducted  by  Rocketdyne.   Among 
the  various  aspects  of  large  excimer  lasers  considered  in  the 
Rocketdyne  reports  was  beam  quality  as  influenced  by  non- 
uniformity  of  index  of  refraction.   The  adverse  impact  of 
poor  beam  quality  upon  the  overall  system  was  stressed. 

For  military  applications  the  optimum  pulse  duration 
should  be  a  few  microseconds  [13-15].   The  pulse  duration 
for  some  military  applications  is  in  contrast  to  the  desired 
pulse  duration  of  a  few  nacoseconds  for  inertial  confinement 
to  achieve  thermonuclear  fusion  [1].   Even  if  a  pulse  dura- 
tion of  10  microseconds  is  used,  the  motion  of  the  waves, 
e.g.,  cathode  shock  wave  and  anode  wave,  in  the  gas  will  be 
very  small,  i.e.,  of  the  order  of 

I     =      a  t   =  (300  m/sec)  (10~5  sec)  =  3  mm    (1) 

where  a  is  the  acoustic  wave  speed  and  t  is  the  pulse  dura- 
tion.  However,  this  fact  does  not  suggest  that  a  pulsed 
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excimer  laser  will  be  free  of  inhomogeneity  problems ;  exten- 
sive experience  with  pulsed  electrical  lasers  operating  in 
the  IR  indicates  that  beam  quality  may  be  a  problem.   To 
achieve  high  average  power,  a  large  pulse  repetition  fre- 
quency, PRF ,  must  be  used  which  requires  a  flowing  laser. 
The  higher  the  PRF,  the  faster  the  gas  flow  is. 

Gas  flow  in  a  high  power,  high  PRF,  excimer  laser  will 
probably  be  subsonic.   In  subsonic  flow  [16-13],  disturbances 
can  propagate  upstream.   Even  if  the  upstream  base  flow  is 
perfect,  the  gas  moving  out  of  the  cavity  after  a  pulse  can 
cause  An/n  within  the  cavity;  n  is  index  of  refraction,  and 
An  is  spatial  variation  of  index  of  refraction.   A  large 
flush  factor  eases  the  problem  of  a  nonuniform  base  flow; 
however,  high  efficiency  demands  a  low  flush  factor. 

Nonuniformity  of  the  lasing  medium  may  have  adverse 
effects  distinct  from  the  deterioration  of  beam  quality. 
The  work  of  Nerheim  and  Chen  [19]  at  JPL  indicates  different 
levels  of  performance  depending  on  homogeneity  of  the  gas . 
Pulses  into  a  homogeneous  medium  resulted  in  one  level  of 
power  output;  pulses  into  an  inhomogeneous  medium  yielded 
another  level  (lower  naturally)  of  performance.   Details  of 
the  physics  involved  in  this  observation  were  not  unraveled. 
The  work  by  Nerheim  and  Chen  was  with  an  atomic  copper  laser. 
Consequently  the  results  are  not  directly  transferable  to 
rare  gas-halide  lasers. 

Despite  the  many  unsolved  problems  regarding  excimer 
lasers,  its  use  as  a  potential  weapon  constitutes  a 
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revolution  greater  than  the  leap  from  naval  guns  to  guided 
missiles.   Some  of  the  conceivable  missions  for  the  excimer 
laser  are  ballistic  missile  defense  (both  ground-based  silo 
defense  and  space-based  destruction  of  ICBMs  in  the  boost 
phase),  space-based  antisatellite  defense,  and  antiship 
missile  defense  (ASMD)  [20].   The  latter  has  the  most  pro- 
mise of  relatively  near-term  development.  Also,  the  radia- 
tion from  certain  excimer  lasers  has  good  propagation  in 
sea  water. 

Briefly,  the  laser  offers  the  advantage  of  depositing 
an  extremely  high  energy  on  a  small  area  of  a  target.   This 
destructive  energy  is  delivered  at  the  speed  of  light  (300 
million  meters  per  second  which  is  ideal  for  a  defensive 
system  against  missiles);  whereas,  a  conventional  missile 
travels  to  its  target  at  a  speed  of  a  few  times  the  speed 
of  sound  (roughly  1500  meters  per  second) . 

As  alluded  to  earlier,  beam  quality  has  an  important 
impact  upon  the  overall  laser  system.   Improvements  in  beam 
quality  can  be  used  to  decrease  overall  laser  size,  shorten 
time  for  melting  a  metal,  or  decrease  laser  run  time  to 
perform  a  certain  task.   A  well  know  fact  is  that  beam 
quality  depends  on  favorable  illumination  of  the  exit 
aperture;  favorable  illumination  implies  near  constant 
phase  and,  to  a  lesser  degree,  near  constant  amplitude. 
Variation  of  the  index  of  refraction  within  the  cavity 
causes  spatial  variation  in  phase  at  the  laser  exit.   This 
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nonuniform  phase  causes  the  peak  axis  intensity,  Iff,  to 
decrease  according  to 


J£   =  exp  [-(2T^rms/A)2]  (2) 

o 


where  X  is  wavelength.   The  rms  phase,  4)rms »  is  obtained  by 
squaring  and  averaging  the  phase  shifts  across  the  aperture. 
Equation  (2)  is  the  Strehle  Ratio  [21].   The  spatial  phase 
relation  is 

<j>  »  /  nds  (3) 

where  ds  is  an  incremental  distance  along  a  ray  in  the 
laser  cavity.   The  initial  efforts  to  limit  phase  distortion 
considered  the  index  of  refraction  in  terms  of  the  Gladstone- 
Dale  formulation 

.  n  =  1  +  Gp  (4) 

where  p  is  gas  density,  and  G  is  the  Gladstone-Dale  constant; 
hence  there  is  concern  for  flow  and  acoustic  effects.   Hogge 
and  Crow  [22]  state  that  flow  uniformity  and  reduced  acoustic 
effects  have  achieved  Ap/p  =  5x10   ,  which  is  the  required 
flow  uniformity. 

It  has  been  suggested,  however,  that  the  requirement  for 
flow  uniformity  should  be  stated  in  terms  of  An/n  due  to 
nonuniform  pumping  in  the  cavity  [23].   In  a  typical  flow 
mixture,  the  constituents  of  the  medium  (Ne:Xe:NF~)  are 
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present  in  a  fixed  ratio  (750:3:1)  [24].   Even  in  a  gas  with 
fixed  chemical  composition,  pressure  and  translational  tem- 
perature, and  hence,  fixed  Ap/p,  An/n  can  vary  according  to 
the  ionic  and  electronic  excitation  of  the  medium  caused  by- 
pump  ing . 

Consider  the  Ladenburg  formulation  of  the  Kramers -Kronig 
dispersion  relation  [25] 


2M   x-^   N.  A,2.   X2  N,g. 


51      j   k        kj 


(5) 


e 

electron  charge 

c 

speed  of  light 

m 
e 

mass  of  electron 

N 

total  number 

A,  .  wavelength  for  transition  between  k  and  j 

A  wavelength  of  interest  at  which  n  is  derived 

f,  .  is  oscillator  strength  for  transition  from  k  to  j 

g.  is  degeneracy  of  j    energy  level 

t*V» 

N.    population  in  j    energy  level 


There  are  I   distinct  species  indicating  specific  electronic 
or  ionic  excitation  present  in  the  cavity.   Index  j  sums 
over  lower  levels  and  k  over  upper  levels  of  possible  con- 
figurations .   Identify  the  mole  fraction  as 
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X,  =  gi  (6) 


the  fraction  of  the  total  molecular  population  represented 
by  species  I.      Finally,  define  a  quantity 


2   2 
2     ^       xf.X  N.  g. 

e     j    k.         tcj 


In  terms  of  the  quantities  defined  above 

n-l  -  H   £  X,  T,  (8) 

I 

or  the  index  of  refraction  in  the  cavity  is  the  sum  of  the 
contributions  of  all  the  species  present  in  the  cavity. 
Since  U     =  NX., 

n-l  =  N1T1  +  M2T2  +  .  .  .  (9) 


The  partial  density  is  given  by  p   =  N  m  where  m   is  the 
molecular  weight  of  I        species .   In  terms  of  partial 
density,  equation  (9)  becomes 


n-l  =  p1T1/m1  +  p2T0/m2  +  ...         (10) 


Introduce  the  mass  fraction  for  the  gas  mixture  Y.  =  P9/p. 
where  p  is  the  density  of  the  gas  mixture.   The  result  is 

n-l  =  p[Y1T1/m1  +  Y2T2/m2  +  ...]      (11) 
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Equation  (11)  shows  that  n-1  is  proportional  to  p;  however, 
equation  (11)  further  demonstrates  that  n-1  depends  on  mass 
fraction  and  index  of  refraction  of  each  species .   Even  with 
uniform  flow,  the  spatial  phase  relation  across  the  cavity 
will  not  be  uniform.   Accurate  modeling  of  the  phase  effects 
requires  knowledge  of  the  species  populations  and  the  con- 
tribution to  the  index  of  refraction  of  the  species  in  the 
cavity.   This  thesis  will  explore  a  model  for  determining 
the  species  populations  in  a  laser  cavity  and  apply  a  model 
originally  developed  at  the  Naval  Research  Laboratory  (NRL) 
to  calculate  the  populations  of  the  different  configurations 
of  excited  Neon  present  in  a  XeF  laser. 

First,  in  section  II,  a  review  of  electronic  transition 
gas  laser  kinetics  is  discussed.   Second,  in  section  III,  a 
simple  model  for  an  excimer  (KrF)  laser  is  developed.   Third, 
in  section  IV,  the  computer  program  developed  and  written  at 
NRL  by  Johnson,  Palumbo ,  and  Hunter  [26]  is  discussed. 
Fourth,  in  section  V,  a  formulation  of  input  reactions  to 
the  NRL  program  to  calculate  the  populations  of  excited 
neon  configurations  is  developed.   Fifth,  in  section  VI, 
calculation  of  the  index  of  refraction  is  performed.   Sixth, 
in  section  VII,  the  calculation  of  phase  shift  is  performed. 
Lastly,  in  section  VIII,  the  results  are  discussed. 
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II.   KINETICS 

Electronic  transition  gas  lasers  (i.e.  exciplex  or 
excimer)  generally  involve  complex,  branching  chains  of 
energy  flow  from  initial  pump  excitation  of  background  gas 
to  eventual  extraction  of  laser  radiation.   Figure  2  illus- 
trates major  energy  flow  pathways  in  e-beam  pumped 
Ne/Xe/NF-,  mixtures  .   These  energy  chains  include  processes 
such  as:   energy  deposition  by  electron  ionization  and 
excitation;  energy  transfer  by  various  collisional 
mechanisms,  including  neutralization,  displacement  and 
"harpooning",  leading  to  formation  of  the  upper  laser 
level;  and  collisional  mechanisms  which  intercept  excitation 
energy  before  it  can  reach  the  upper  laser  level,  or  which 
quench  that  level  before  stimulated  emission  occurs.   Also, 
one  must  consider  extraction  of  stimulated  emission  includ- 
ing absorptions  by  various  processes  in  the  gas  and  reso- 
nant reabsorption  by  the  lower  laser  level  (if  there  is 
one) ,  laser  cavity  oscillation,  and  coupling  through  output 
mirrors.   A  time-dependent  model  of  the  kinetics  of  such  a 
system  must  define  and  update  species  densities  for  a  large 
number  —  typically  between  20  and  69--gas  components,  includ- 
ing ground  states,  several  excited  states  and  ions  of  atoms 
and  molecules,  plus  electrons  and  photons. 

The  time  evolution  of  these  components  is  followed  by 
the  solution  of  a  set  of  coupled  first  order  ordinary 
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Figure  2.   Major  energy  flow  pathways  in  e-beam  pumped 

Ne/Xe/NF3  mixtures  (Reproduced  from  Huestis  [27]) 

Note:       Energy  input  from  the  electrons  comes  in  along 
the  pump  arrows.   The  species  in  each  box,  plus 
the  component  along  the  reaction  arrow,  yield 
the  result  in  the  next  box.   The  wavelengths  of 
the  emitting  species  are  indicated.   They  may  be 
subject  to  quenching  by  NF3  or  Xe  (not  shown). 
NeF*  is  presumed  to  predissociate . 
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differential  equations;  one  equation  is  required  for  each 
gas  component.   To  the  set  of  equations  defining  species 
densities  are  added  further  equations  describing  the  evolu- 
tion of  the  photon  field,  the  gas  temperature,  and  (in  some 
models)  coupled  electrical  circuit  parameters.   The  species 
equations  can  be  generally  described  as  equating  the  rate 
of  change  of  species  density  to  a  sum  of  contributing  for- 
mation terms  less  a  sum  of  depletion  terms  as 


dNi    V-         V- 


Here  N.  is  a  particular  species  density,  the  F . .  are  the 

formation  terms,  and  the  D.,  the  depletion  terms.   For 

lk 

example,  consider  a  simple  chemical  reaction  of  the  form 
A  +  B  -*-  C  +  D  whose  forward  rate  constant  is  k;  if  the 

species  densities  are  expressed  in  cm   ,  the  units  of  k  are 

3 
cm  /sec.   For  this  reaction,  the  rate  equations  defining 

species  C  and  D  would  contain  the  formation  term 

F  «  k[A][B]  ;  (13) 

at  the  same  time,  the  equations  defining  species  A  and  B 
would  contain  the  depletion  term 

D  -  k[A][B]  .  (14) 

Recall  that  square  brackets  for  a  chemical  symbols  means 

3 
number  density  of  the  chemical  with  dimensions  1/cm  . 
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Consider  equation  (12)  again.   Species  i  is  formed  in 
j  different  reactions.   Likewise,  species  i  is  depleted  in 
k  different  reactions . 
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III.   A  SIMPLE  KINETICS  MODEL  FOR  KrF 

In  a  thorough  model  of  an  excimer  laser  a  large  number 
of  kinetic  reactions  are  necessary;  see,  for  instance,  the 
185  reactions  listed  in  the  output  of  the  NRL  program  for 
XeF  (Appendix  A) .   The  reactions  are  used  to  obtain  the  set 
of  coupled  first  order  differential  equations .   The  solution 
to  this  set  of  equations  is  very  involved  and  tedious 
because  of  the  difficulty  in  handling  the  tangle  of  reactions 
Therefore,  only  a  few  of  the  kinetic  reactions  in  a  KrF 
laser  are  considered.   The  object  of  the  simple  computer 
model  was  to  obtain  a  qualitative  insight  to  excimer  laser 
kinetics . 

The  model  uses  a  He/Kr/F2  mixture  (93.1/1.77/0.13)  with 
He  as  diluent.   Figure  3  illustrates  a  simplified  e-beam 
pumped  KrF  laser  geometry. 


ELYTRON     BEAM 

I   I   I 


>^   wave  path 


Point  of  Stimila-ted 

Emission 

LASER  CAVITY 


<-CAIN  LENGTH  (L) 


H 


Figure  3.   Example  of  e-beam  pumped  KrF  laser  geometry 
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The  first  step  in  the  kinetics  chain  leading  to  extrac- 
tion of  laser  radiation  is  the  deposition  of  energy  in  the 
diluent  gas  by  either  photons  or  electrons.   In  this  model 
electrons  are  used.   The  electrons  '  ionize  the  Kr  and 
dissociate  the  F« 

kl   + 
e'  +  Kr   +  Kr  +  e  +  e'  (15) 

where  e1  is  a  primary  electron,  e  is  a  secondary  electron 
resulting  from  ionization  of  Kr,  and  k,  is  the  forward 
rate  constant.   The  secondary  electron  has  sufficient 
energy  to  dissociate  and  electronegatively  attach  to 
atomic  F . 

k2 

e  +  F2   -   F  +  F  (16) 

In  addition  dissociation  of  F~  by  the  primary  electrons 
occurs . 

k3 

e'  +  F2   ->  F  +  F~  (17) 

Having  produced  the  ion  pair  states,  recombination  of 
positive  rare  gas  ions  and  readily  formed  halide  ions  yields 
the  excited  state  (KrF*) . 

+      - 

Kr  +  F   +  M   ■*   KrF*  +  M  (18) 

The  KrF*  can  produce  a  photon  by  stimulated  emission. 

k5 

KrF*  +  hv   ■+   KrF  +  2hv  (19) 
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Also,    recombination   of  F   can   occur 

k6 

F   +  F+M      ■*•        F2+M  (20) 

Fluorescence  decreases  KrF*  population  by 


T 

KrF*   ■+  KrF  +  hv  (21) 


where  t  is  the  radiative  lifetime. 

Reactions  (15)  to  (21)  are  the  only  reactions  considered 
in  this  model.   None  of  the  quenching  processes  of  KrF* 
have  been  considered,  such  as 


KrF*  +  2Kr   ■+  Kr£F*  +  Kr  (22) 


nor  have  many  other  reactions  [2]  and  [28]. 

The  reactions  considered  in  this  model  are  used  to  find 
the  formation  terms  equation  (13)  and  depletion  terms 
equation  (14).   For  example,  examining  reaction  (15)  the 
formation  term  for  the  concentration  of  e  and  Kr   is 
determined . 

Formation   e  =  k-Je'HKr]  (23) 

Formation  Kr+  =  k1[e'][Kr]  (24) 

3 
where  [  ]  denote  concentration  in  specie  number/cm  . 

The  depletion  term  for  concentration  of  e  from  reaction  (16) 

is 
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Depletion   e  =  k2[e][F2]  (25) 

and  for  Kr   from  reaction  (18) 

Depletion  Kr+  =  k4 [Kr+] [F~] [M]  (26) 

Applying  equation  (12)  ,  that  is  ,  equating  the  rate  of  change 
of  e  concentration  to  the  sum  of  contributing  formation 
terms  equation  (23)  less  the  sum  of  depletion  terms  equation 
(25)  yields 

|^i    =  k^e'HKr]  -  k2[e][F2]  (27) 

Similarly,  combining  equations  (24)  and  (26) 

d[Kr  ]   =  k^e'HKr]  -  k4[Kr+][F"][M]  (28) 

Continuing  the  process 

^J-    =  k2[e][F2]  +  k3[e'][F2]  -  kg[F][M]     (29) 

1^1    =  k2[e][F2]  +  k3[e'][F2]  -  k4 [Kr+] [F~ ] [M] 

(30) 

d[KrF*]  =  k4[Kr+][F"][M]  -  k5[hv][KrF*]  -  [KrF*]/x 

(31) 
d[KrF] 


dt 


-  k5[KrF*][hv]  +  [KrF*]/i  (32) 


are  obtained. 

In  reaction  (15)  the  e'  is  the  e-beam  input  and  there- 
fore is  a  function  of  time.   The  model  assumes  the 
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3 

concentration  of  e'  is  initially  1.9E13  electrons /cm   then 

becomes  zero  after  about  177  nanoseconds.   The  secondary- 
electron  mean  energy  is  10  eV . 

The  rate  constants  (k)  for  the  reactions  were  taken 
from  scientific  articles  [2,  29,  and  30]  or  derived  as  in 
the  case  of  k,-  .   The  values  of  the  rate  constants  are  listed 
in  Table  II.   The  radiative  lifetime,  t,  is  equal  to  15 
nanoseconds  from  reference  [2]. 


Table  II.   Rate  Constants 

Rate  Constant Value Reference 

k1            5.0E-09  cm3/sec  29 

k2            2.0E-09  cm3/sec  30 

k3            5 .0E-10  cm3/sec  30 

k4            1.0E-25  cm6/sec  2 

t,                                i    ctt    nc         3,  Derived   in 

k,            1.5E-05  cm  /sec  .     , .   „ 

5  Appendix  3 

k6            1.0E-31  cm6/sec  2 


Equations  (27)  to  (32)  can  be  solved  to  calculate  the 
concentrations  as  a  function  of  time.   Multiplying  these 
equations  by  dt  yields  a  form  for  numerical  solution  on  the 
HP9830.   Equations  (27)  to  (32)  become 

d[e]  =  (k^e'HKr]  -  k2[e]  [F2]  )dt  (33) 

d[Kr]  =  (k^e'HKr]  -  k4  [Kr+]  [F~  ]  [M]  )  dt  (34) 
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[e] 

=   0 

[F] 

=   0 

[F"] 

=   0 

[Kr+] 

=   0 

d[F]  =    Ck2[e][F2]    +k3[e'][F2]    -k6[F][M])dt  (35) 

d[F"]         =    (k2[e][F£]    4-k3[e'][F2]    -    k4[Kr+] [F" ] [M] ) dt 

(36) 

d[KrF*]  -  (k4[Kr+] [F~] [M]  -  k5[hv][KrF*]  -  [KrF*]/x)dt 

(37) 

d[KrF]   =  (k5[KrF*][hv]  +  [KrF*]/x)dt  (38) 

respectively . 

The  initial  conditions  for  the  set  of  equations  are: 

(39) 
(40) 
(41) 
(42) 

[KrF]   =  0  (43) 

[KrF*J  =  0  (44) 

Several  constants  are  specified  as  follows: 

[hv]    =  4.0E11  (calculated  in  Appendix  C)  (45) 

[F2]    -  4.0E16)  (46) 

/  Assumed  values  for 
[Kr]    =  5.4E17\the  He/Kr/F2  mixture  (47) 

(  (93.1/1.77/0.13)  at 
[He]    =  3.0E191  STP .  (48) 

Initial  T  =  0 

Final  T   =  250  nanoseconds 

Delta  T   =  (Final  T  -  Initial  T)/256 

A  computer  program  was  written  using  equations  (33)  to 
(38)  and  the  initial  conditions  and  appropriate  constants 
given  by  equations  (39)  to  (48).   Figures  4  and  5,  which  are 
plots  obtained  from  the  HP9330  program,  are  the  concentrations 
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3 
in  specie  number/cm  versus  time  from  0  to  250  nanoseconds. 

The  plots  in  Figures  4  and  5(b)  indicate  that  some  species 
densities  rise  extremely  steeply  in  the  first  tenth  of  the 
simulation  time;  later  in  the  pulse  (after  about  50  nano- 
seconds) species  densities  tend  to  flatten  out.   See  Figure  4 
Once  the  e-beam  terminates ,  the  decay  time  is  approximately 
50  nanoseconds.   In  Figure  5(a)  and  (b)  other  species 
densities  increase  gradually  and  flatten  out.   The  results 
of  this  model  are  in  qualitative  agreement  with  the  NRL 
simulation  contained  in  reference  [26].   The  NRL  simulation 
limited  the  change  in  concentration  to  5  percent  per  time- 
step  to  obtain  sufficient  accuracy  [26].   The  question  of 
timesteps  for  the  simple  model  was  avoided;  sufficiently 
small  steps  were  used  throughout  the  program. 

Table  III  lists  the  variables,  symbols,  and  equations 
used  in  HP9830  computer  program. 

Table  IV  lists  the  KrF  Computer  Program  and  Table  V 
lists  the  output. 
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Figure  4.   Concentration  of  secondary  electrons  (e) ,  Kr '  and 
KrF,v ,  molecules / cm^  ,  as  a  function  of  time 
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Figure  5 


Concentration  of  F,  F  ,  and  KrF  molecules/cm  , 
as  a  function  of  time. 
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Table  III.   Variables,  Symbols ,  and  Equations  Used  in 
HP9830  Computer  Program  for  KrF  Laser. 


Variables  and  Symbols  Used  in  the  Program 
Text    Meaning  or  Value  (for  constants)     Units 
primary  electron 
secondary  electron 
Fluorine  atom 
Fluorine  molecule 

Fluorine  atom  with  extra  electron 
Krypton 

Krypton  minus  an  electron 
Third  body 
Krypton  Fluoride 
Krypton  Fluoride  excited 

photon  joules 

3 
Rate  constant  cm  /sec 

3 
Rate  constant  cm  /sec 

3 
Rate  constant  cm  /sec 

Rate  constant  cm  /sec 

3 
Rate  constant  cm  /sec 

Rate  constant  cm  /sec 

Radiative  lifetime  sec 

Index;  1=1  for  initial  conditions 

Time  increment  sec 


El 

e' 

E 

e 

F 

F 

G 

F2 

H 

F" 

K 

Kr 

L 

Kr+ 

M 

M 

X 

KrF 

Y 

KrF* 

P 

hv 

Kl 

kl 

K2 

k2 

K3 

k3 

K4 

k4 

K5 

k5 

K6 

k6 

Tl 

T 

I 

- 

T 

t 
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Table  III  (Continued).   Variables,  Symbols,  and  Equations 

Used  in  HP9330  Computer  Program  for  KrF  Laser. 

Equations  Used  in  the  Program 
Line  Number  Equation  Number 

860  33 

370  34 

880  35 

890  36 

900  45 

910  37 

920  38 
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IV.   XeF  LASER  COMPUTER  PROGRAM  FROM  URL 

In  order  to  calculate  the  excited  states  of  neon  in  a  XeF 
laser  a  complex  kinetics  program  is  required.   A  kinetics 
scheme  for  the  XeF  laser  developed  at  the  Naval  Research  Lab 
was  used.   It  was  written  by  Dr.  Louis  J.  Palumbo  [26  and  31] 
who  assisted  greatly  in  modifying  it  for  use  on  the  IBM  360/67. 

The  computer  model  consists  of  two  parts:   first,  an  e- 
beam  deposition  section,  in  which  the  kinetic  energy  of  near 
relativistic  electrons  is  used  to  ionize  and  excite  the  laser 
gas,  and,  second,  a  chemical  kinetics  section,  in  which  the 
energy  deposited  in  the  gas  cascades  down  by  various  processes 
to  the  upper  laser  level  and  is  extracted  as  laser  radiation. 

One  input  to  the  NRL  computer  model  is  a  list  of  reactions, 
e.g.  in  the  format  of  equations  (16)  to  (21).   Also  an  input 
is  a  list  of  reaction  rate  coefficients.   The  computer  program 
sorts  the  reactions  and  forms  formation  terms,  F,  and 
depletion  terras,  D.   Formation  and  depletion  terms  were 
discussed  earlier  in  connection  with  equations  (13)  and  (14) . 
Continuing,  the  computer  program  selects  the  appropriate  F 
and  D  terms  to  form  the  rate  equation  of  equation  (12).   Next, 
the  computer  program  multiplies  the  rate  equation  by  a  variable 
times tep.   The  same  approach  was  used  in  the  simple  computer 
program;  see  equations  (33)  to  (38) . 
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An  advantage  of  the  NRL  program  is  the  ability  to  change 
the  number  of  reactions.   The  computer  program  does  not  change; 
only  the  input  changes.   The  differential  equations  are 
generated  internal  to  the  program.   This  method  makes  it  much 
easier  to  alter  the  set  of  equations  and  greatly  reduces  errors 
attendant  to  such  alterations  which  arise  with  ease  in  the 
complex  tangle  of  reactions.   The  detailed  explanation  of  the 
computer  model  is  contained  in  references  [26]  and  [31]. 
Appendix  A  gives  a  list  of  reactions  used  to  generate  the 
population  distribution  within  neon  energy  levels. 


V.   CALCULATION  OF  NEON  EXCITED 
STATE  POPULATIONS 

The  excited  states  of  neon  are  represented  in  Figure  6. 
The  electronic  configurations  of  interest  to  this  study  are 
the  3s,  3p,  4s,  3d,  4p ,  5s,  4d,  4f,  5p,  5d,  and  6s. 

In  order  to  calculate  the  populations  of  these  particular 
configurations,  one  needs  to  know  the  reactions  involved  and 
the  rate  constants.   For  example, 

ki 

Ne  +  e-beam  -*>  Ne(3s) 

k2 

Ne  +  e-beam  ■*  Ne(3p) 


kll 
Ne  +  e-beam  -  Ne(6s) 


These  reactions  could  then  be  input  to  the  NRL  XeF  pro- 
gram to  calculate  the  populations  of  the  excited  configura- 
tions of  neon.   After,  an  extensive  literature  search  (Optics 
and  Spectroscopy,  The  Journal  of  Chemical  Physics,  Applied 
Physics  Letters,  Physical  Review  A,  Journal  of  the  Optical 
Society  of  America,  Chemical  Abstracts,  IEEE  Journal  of 
Quantum  Electronics,  etc)  rate  constants  for  the  individual 
excited  configurations  i.e.  3s,  3p ,  ...  were  not  found. 

Therefore,  it  was  necessary  to  make  a  number  of  assump- 
tions.  The  procedure  followed  was  to  extract  all  reactions 
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Figure    6.      Energy-Level   diagram  of  Neon    (Reproduced   from 
Pressley    [ 32] . 
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L.20E-13 

(49) 

0.925E-18 

(50) 

0.275E-18 

(51) 

from  the  NRL  program  involving  the  neon  excited  state  (Ne*) . 
Then,  these  reactions  were  subdivided  into  two  groups,  the 
Ne(3s)  and  all  higher  configurations,  which  are  denoted  by 
NE*1  and  NE*Z  respectively.   The  population  distribution 
within  NE*Z  was  calculated  assuming  a  Boltzmann  distribution. 
The  list  of  reactions  which  were  used  as  modified  input 
to  the  NRL  program  are  listed  in  Table  VI  in  the  following 
format : 

Rate  Constant 

Ne  +  e'  *  Ne*  +  e* 

Ne  +  e'  ■+  Ne*l  +  e' 
Ne  +  e'  t-  Ne*Z  +  e' 

To  compare  equations  (49)  to  (51)  with  the  first  three 
reactions  of  Table  VI  note  that  e'  is  EBEM  in  FORTRAN  nota- 
tion.  The  other  symbols  are  obvious.   The  first  reaction 

yielding  Ne*  is  divided  into  the  Ne*l  and  Ne*Z  reactions. 

o 
A  computer  run  with  the  current  density  equal  to  10  A/ cm 

and  a  computer  run  with  the  current  density  equal  to  11 

2 

A/ cm  calculated  the  populations  for  Ne*l  (3s)  and  Ne*Z 

(all  levels  greater  than  3s). 

The  rate  constants  for  the  reactions  were  determined  by 
taking  a  fraction  of  the  total  rate  constant.   For  example, 
in  Table  VI,  for  the  second  set  of  reactions,  the  rate 
constant  for  NE*1  is  0,8  and  NE*Z  is  0.2  of  the  rate  con- 
stant for  NE*.   These  fractions  of  the  rate  constant  for 
NE*  are  noted  in  parentheses  following  the  reactions.   The 
choices  of  the  fractional  values  were  arbitrary. 
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TABLE  VI.   Input  Reactions  to  the  NRL  Program. 


* 


Reaction 


Ps.ate   Constant 


NE 

+ 

EBEM 

■> 

NE* 

+ 

E3EM 

1.20 

E-18 

NE 
NE 

+ 
+ 

EBEM 
SBEK 

NE*1 
NE*Z 

+ 
+ 

EBEM 
E3EM 

0.925E-18 
0.275E-18 

NE* 

+ 

77 

-V 

NE 

+ 

E- 

2.00 

E-10 

NE*1 

NE*Z 

+ 
+ 

E- 

NE 
NE 

+ 
+ 

E- 
E- 

(.3) 
(.2) 

1.60 
0.40 

E-10 

E-10 

NE2+ 

+ 

E- 

-»■ 

NE* 

+ 

NE 

3.50 

E-08 

NE2+ 
NE2+ 

+ 
+ 

E- 
E- 

NE*1 
NE*Z 

+ 
+ 

NE 

NE 

(.9) 
CD 

3.15 
0.35 

E-08 
E-08 

NE3+ 

+ 

E- 

->- 

NE* 

+ 

NE 

+ 

NE 

3.50 

E-08 

NE3+ 
NE3+ 

+ 
+ 

lit- 

NE*1 
NS*Z 

+ 
+ 

NE 

NE 

+ 
+ 

NE 
NE 

(.95) 
(.05) 

3.325E-08 
0.175E-08 

XENE+ 

+ 

E- 

-»■ 

NE* 

+ 

XE 

6.00 

E-08 

XENE+ 
XENE+ 

+ 
+ 

E- 
E- 

NE*1 
NE*Z 

+ 
+ 

XE 
XE 

(.6) 
(.4) 

3.60 
2.40 

E-08 
E-08 

NE* 

+ 

NE 

+  NE 

-► 

NE2* 

+ 

NE 

5.0 

E-34 

NE*1 
NE*Z 

+ 
+ 

NE 
NE 

+  NE 
+  NE 

->- 

NE2* 
NE2* 

+ 
+ 

NE 
NE 

(.8) 
(.2) 

4.0 
1.0 

E-34 
E-34 

NE* 

+ 

XE 

-V 

XE+ 

+ 

E- 

+ 

NE 

7.40 

E-ll 

NE*1 
NE*Z 

+ 
+ 

XE 
XE 

->- 
-> 

XE+ 
XE+ 

+ 
+ 

E- 
E- 

+ 
+ 

NE 
NE 

(.6) 
(.4) 

4.44 
2.96 

E-ll 
E-ll 

NE* 

+ 

XE 

->- 

NEXE+ 

+ 

E- 

2.20 

E-ll 

NE*1 
NE*Z 

+ 
+ 

XE 
XE 

-> 
->■ 

NEXE+ 
NEXE+ 

+ 
+ 

E- 

E- 

(.6) 
(.4) 

1.32 
0.88 

E-ll 
E-ll 
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TABLE  VI  (Continued) .   Input  Reactions  to  the  NRL  Program. 


Reaction 

Rate  Constant 

NE* 

+  F2 

-*  NEF* 

+  F 

4.10  E-10 

NE*1 
NE*Z 

+  F2 
+  F2 

-►  NEF* 
■>  NEF* 

+  F 
+  F 

(.8) 
(.2) 

3.28  S-10 
0.62  E-10 

NE* 

+  NF3 

■+   NEF* 

+  NF2 

1.05  E-10 

NE*1 
NE*Z 

+  NF3 
+  NF3 

+   NEF* 
+   NEF* 

+  NF2 
+  NF2 

(.8) 
(.2) 

0.84  S-10 
0.21  E-10 

NE* 


+  NF2 


NEF* 


+  NF 


1.0   E-10 


NE*1 

+  NF2 

•*■  NSF* 

+  NF 

NE*Z 

+  NF2 

H.  NEF* 

+  NF 

(.3) 

0.8 

S-10 

(.2) 

0.2 

E-10 

NE*    +  NF 


NEF* 


+  N 


1.0   E-10 


NE*1 

+  NF 

+  NEF* 

+  N 

NE*Z 

+  NF 

+   NEF* 

+  N 

(.8) 

0o8 

E-10 

(.2) 

0.2 

E-10 

NE2*   +  HN352L 


NE* 


+  NE   +  ABSL 


1.0   E-18P 


NE2*  +  HN35  2L 
NE2*  +  HN352L 


NE* 

+  NE 

+  ABSL 

(.9) 

0.9 

E-18P 

NE* 

+  NE 

+  ABSL 

(.1) 

0.1 

E-18P 

Use  of  capital  letters  more  closely  represents  computer 
printing. 
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However,  the  first  reactions  listed  in  Table  VI  do  not 
have  fractional  values  in  parentheses,  as  the  reaction  rate 
constants  were  determined  using  reference  [9].   The  plot  of 
"Excitation  of  Neon  by  Electrons"  from  page  467  of  refer- 
ence [9]  was  extrapolated  to  10  keV,  and  the  value  for  the 
2p,  state  was  assumed  to  represent  all  excited  configura- 
tions above  the  3s  (NE*1) .   The  difference  between  the  rate 
constant  for  NE*  and  the  2p-,  value  yielded  the  rate  constant 
for  NE*1. 
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VI.   CALCULATION  OF  THE  INDEX  OF  REFRACTION 

The  values  (n-l)/N  where  n  is  the  index  of  refraction 
and  N  is  the  species  population,  have  been  calculated  Refs. 
[33],  [34]  and  [35].   Figure  7  lists  the  values  of  (n-l)/N. 
References  [34]  and  [35]  contain  the  details  of  calculating 
these  values. 

The  index  of  refraction  (n-1)  is  equal  to 

n-l  .  (Igll)  N  +  (^y)  N(3s)  +  (^)  N(Z)   (52) 

ground  n  * 

state 

Using  the  populations  calculated  in  the  NRL  model  for  N, 
N(3s) ,  and  N(Z)  equation  (52)  can  be  solved.   However,  the 
populations  of  individual  electronic  configurations  must  be 
determined. 

One  feature  of  excimer  lasers,  in  contrast  to  the 
vibrational  transition  lasers  of  the  infrared  region,  is 
the  extent  of  randomness  of  populations.   In  the  C02  gas 
dynamic  laser,  mode  temperatures  can  be  assigned  to  the 
various  vibrational  modes .   Due  to  the  short  span  of  excita- 
tion and  power  extraction  in  an  excimer  laser,  the  various 
modes  do  not  have  an  opportunity  to  equilibrate.   This  fact 
is  the  reason  extensive  computer  programs  are  necessary  for 
kinetics.   An  assumption  of  a  Boltzmann  distribution  is 
unfoundedo   However,  in  view  of  the  sparsity  of  reaction 
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rate  coefficients  for  excitation  to  specific  configurations 
a  Boltzmann  distribution  is  assumed.   Assuming  a  Boltzmann 
relation  for  the  upper  levels  the  populations  of  the  3p , 
4s,  3d,  ...  6s  can  be  calculated. 

N(Z)  =  N(3p)  +  N(4s)  +  N(3d)  +  ...  +  N(6s)        (53) 


-(E,  -E~  ) 
4s   3p7 


N(Z)  =  N(3p)  [1  +  —  exp 


'3p 


-<E3d-E3p> 


a3d        kT       -,   ,c,. 
+  r —  exp  . . . ]   (54) 

S3P 


where  g  denotes  the  degeneracy  of  the  configuration  which  is 
equal  to  the  sum  of  the  2J  +  1  terms.   E  is  the  energy  sub- 
scripted by  the  configuration » 
Solving  for  N(3p) 

N(Z) 

N(3P)  =  _/E   .tt   S (55) 

CE4s  -y 

1  +  ^exp    kT 
g3p 

All  the  quantities  on  the  right  side  of  equation  (55) 
are  known  or  can  be  calculated.   Table  VII  lists  the  values 
of  the  degeneracies  and  energies  for  the  configurations. 
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TABLE  VII.   Degeneracy  and  energy  values  for  the 
configurations . 


Electronic 
Configuration 

K(S  2J 

+1) 

E  (cm-1) 

3s 

12 

134674.6 

3p 

36 

150197.8 

4s 

12 

158950.1 

3d 

73 

161874.4 

4p 

36 

163180.3 

5s 

12 

166123.2 

4d 

73 

166500 

4f 

84 

166562 

5p 

36 

167812.5 

5d 

73 

168000  . 

6s 

12 

168750 

The  temperature  T  is  calculated  as  follows  using  the 
Boltzmann  relation 


-(E,  -E,  ) 
cr         3s   3p7 

N(3s)    s3s        kT  ,_K. 

NT2pT  =  g^  exp  (56) 


Solving  for  T 

*n  rNQs)  §2Pl  =  "(E3s"E2p) 


51 


T   -   2a ^ (57) 

k  *n  [mm  !2?] 


NT25T 


s3s 


Substituting  the  value  of  T  into  equation  (55)  yields  a 
value  for  IJ(3p).   To  calculate  N(4s)  one  uses 


-H      -E   ) 

E  Ws  -3p;j 

N(4s)  =  N(3p)  -=^-  exp      kT  (58) 


>3p 


Similarly  the  populations  of  the  other  configurations 
can  be  calculated  i.e.  IT  (3d)  ,  N(4p)  ,  etc.   Once  the  popula- 
tions of  the  configurations  have  been  found  it  is  a  simple 
calculation  to  obtain  n-1  applying  equation  (52)  as 

n-1      n-1  n-1 


Ne(3p)  '    lie  (4s)  '  '"   He  (5  s) 

are  given  along  the  bottom  of  Figure  7. 

A  computer  program  to  calculate  the  index  of  refraction 
(n-1)  using  equations  (52^,  (55),  (57),  and  (58)  was 
written. 

In  Figure  7  the  3s  -»■  4p  transition  contains  an  R*.   R* 
indicates  that  resonance  may  occur  and  can  have  a  large 
effect  on  the  index  of  refraction.   The  spikes  in  Figure  8 
are  located  at  the  resonant  wavelengths.   Resonance  is 
discussed  in  reference  [36]. 
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The  plots  in  Figure  9  are  of  the  calculated  indices  of 
refraction  (n-1)  versus  time  (nanoseconds).   Figure  9(a) 
is  for  a  laser  wavelength  away  from  3s  ■*  4p  (no  resonance) 

for  both  low  and  high  power.   The  current  density  for  low 

2  2 

power  is  10  A/ cm  ,  and  the  high  power  is  11  A/ cm  .   For  a 

nonresonant  wavelength,  the  change  in  n-1  is  in  the  fourth 

decimal  place. 

The  resonance  effect  was  taken  into  account  in  a  second 
computer  run;  see  Figure  9(b).   For  a  resonant  laser  wave- 
length, the  contribution  to  the  index  of  refraction  due  to 
Ne(3s)  was  calculated  assuming  (n-1) /Ne (3s)  equal  to 
-5.0308;  the  corresponding  value  for  Ne(4p)  was  (n-l)/Ne(4p) 
equal  to  1.667508.   Figure  9(b)  is  for  a  laser  wavelength 
near  the  3s  ■>  4p  resonance  for  low  and  high  power.  The 
resonance  has  a  marked  effect  on  the  index  of  refraction 
in  comparison  to  Figure  9(a)  as  evidenced  by  the  ordinates 
of  the  plots.   Also,  Figure  9  reveals  that  increasing  the 
power  decreases  the  index  of  refraction. 

The  right  hand  ordinate  scale  for  Figure  9(a)  is  the 
ratio  (n-1) /(n-1)  ,  where  (n-1)   is  the  value  at  zero  time. 
The  scale  emphasizes  the  fact  that  the  relative  changes  are 
in  the  fourth  decimal. 

Table  VIII  lists  the  variables,  symbols,  and  equations 
used  in  the  HP9830  program  to  calculate  n-1.   mable  IX  lists 
the  HP9330  program  used  to  calculate  n-1  and  Tables  X  and  XI 
list  sample  output  for  low  and  high  power  respectively. 
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(a)  Index  of  refraction  (n-1)  at  laser  wavelength  away  from  3s  -*•  4p 
resonance  for  high  and  low  power. 
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(B)     Index  of  refraction  (n-1)   at  laser  wavelength  near  3s  ■*•  4p  resonance 
for  high  and  low  power. 


Figure   9.      Calculated   index  of   refraction    (n-1)    as    a    function 
of   time. 
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TABLE  VIII.   Variables,  Symbols,  and  Equations  Used  in 
HP9830  Computer  Program  for  Calculation  of 
n-1. 

Variables  and  Symbols  Used  in  the  Program 

HP9830    Text     Meaning  or  Value  (for  constants)     Units 

degeneracy  (2  2J+1) 

energy  cm 

speed  of  light  cm/sec 

3oltzmann  constant  eV/°K 

Contribution  of  excited  state 

Index  of  the  levels  1  to  12 

Temperature  degrees 

Kelvin 

J       -       Time  increment  index 

Equations  Used  in  the  Program 

Line  Number  Equation  Number 

320  57 

930  56 

1030  54 

1050  55 

1130  53 

1350  52 


G 

CT 
O 

F 

T7 

CI 

C 

K2 

k 

M 

n-l/Ne() 

K 

2p,3s, . . 

K() 

T 

56 


TABLE  IX.   Listing  of  the  HP9330  Index  of  Refraction  (n-1) 


Program. 


10  PRINT 

LOU 

P8WER  10fl/CM+2" 

20  PRINT 

1  0  0 

BATH 

0?  9. 716E+19? 1E- 

50? 1 E - 5  U 

110 

DflTR 

5  0  ?  9  „  7 1 6  E  + 1  9  ?  2 . 

35E+13? 1. 166E+13 

120 

DATA 

100?  9. 716E+19? 5 

i.814E+13»2.882E+13 

130 

DflTR 

150  5  9.  71 SE+  19  j  ;; 

. 709E+13?  4. 034E+13 

140 

DATA 

200?  9. 716E+19?  S 

i  -  8  6  8  E  +  1  3  ?  4 .  8  4  3  E  + 1  3 

150 

DATA 

2  5  0  ?  9 ,  7 1 6  E  + 1 9  ?  ': 

i,  572E+13?5,.  35E+13 

160 

DATA 

300?  9.  716E+19?  ': 

1  „ 998E+13?  5. 66 IE +13 

170 

DATA 

350?  9. 716E+19? 1 

.026E+14?5.351E+13 

ISO 

DATA 

400?  9. 716E+19? 1 

.  0  4  1 E  +  1 4  ?  5  .  9  6  6  E  +  1  3 

190 

DATA 

450?  9, 716E+19? 1 

. 051E+14?  6. 035E+13 

200 

DATA 

500?  9. 716E+19? ] 

. 057E+14* 6. 077E+13 

210 

DATA 

550?  9. 716E+19? 1 

. 06E+14?  6„ 103E+13 

220 

DATA 

600?  9 . 7 1 6  E  + 1 9  ? ] 

.062E+14?  6. 1 18E+13 

2  3  0 

DATA 

6  5  0  ?  9 .  7  1 6  E  + 1 9  ?  ! 

.  U63EH4?  6.  127E+13 

240 

DATA 

700?  9. 716E+19? : 

. 064E+14?  6. 132E+13 

250 

DATA 

750?9.716E+19? 1 

.  ,065E+14?6. 135E+13 

260 

DATA 

3  0  0  ?  9 .  7  1 6  E  +  1 9  ?  : 

. . 065E+14?  6. 137E+13 

270 

DATA 

8  5  0  ?  9  .  7 1 6  E  + 1 9  ?  : 

L.065E+14J6. 138E+13 

2  3  0 

DATA 

900?  9. 716E+19? : 

. . 066E+14?  6. 139E+13 

2  9  0 

DATA 

950?  9. 716E+19? ] 

L . 066E+14?  6. 14E+13 

3  0  0 

DATA 

1 0  0  0  ?  9 .  7 1 6  E  +  1 9 ; 

1  ,,  0  6  6  E  +  l  4  j  6 .  1 4  E  + 1 3 

3  i  0 

DATA 

1050?  9. 716E+19; 

1  „ 866E+14" 6. 14E+13 

3  2  0 

DATA 

1 100?  9. 716E+19' 

1 . 067E+14? 6. 14E+13 

3  3  0 

DATA 

1 150?  9. 716E+19 

1 . 067E+14?  6. 141E+1 

340 

DATA 

1200?  9. 71SE+19 

1.067E+14) 6. 141E+1 

350 

DATA 

1250?  9. 716E+19 

1 . 067E+14?  S. 141E+1 

3  6  0 

DATA 

1300?  9. 716E+19 

1 . 063E+14?  6. 142E+1 

370 

DATA 

1350?  9. 716E+19 

•  1 . 0S8E+14»6. 142E+1 

3  3  0 

DATA 

1400?  9. 716E+19 

< 1 . 069E+14?  6. 143E+1 

3  9  0 

DATA 

1450?9,716E+19 

. 1 . 069E+14?  6. 144E+1 

400 

DATA 

1500?  9. 716E+19 

. i.07E+14>6. 145E+13 

410 

DATA 

1550?  9. 716E+19 

}  1 .  07I-  +  14?  6.  146E+13 

420 

DATA 

1600?  9.  71-6E+19 

•  l„071E-i-14?6„  147E+1 

430 

DATA 

1650?  9. 716E+19 

?  1. 072E+14?6. 143E+1 

440 

DATA 

1700?  9. 716E+19 

j  1 . 073E+14?  6. 149E+1 

450 

DATA 

1750? 9. 716E+19 

< 1 . 074E+14>6. 15E+13 

460 

DATA 

1300?  9.716E+19 

?  1.074E+14j6. 152E+1 

470 

DATA 

1350? 9.71 6E+ 19 

' 1 . 075E+14?  6. 153E+1 

430 

DATA 

1  9  0  0  ?  9  .  7 1 6  E  + 1 9 

? 1 . 077E+14j6. 155E+1 

490 

DATA 

1950?  9,.  716E+19 

? 1.073E+14?  6.  157E+1 

5  0  0 

DATA 

2000?  9. 716E  +  19 

? 1 „ 079E  +  14?  6. 159E+1 

510 

D I M 

TC41  ]?KC41  ]?3C4 

1  Is  ZC  41  1?MC  13]?G[  13 

515 

D I M 

RC41?  12  1?LC4i  ]? 

HC  41  ]?EC  13]?  AC  41?  12 

] 

]?  HE  41?  12] 
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TABL2  IX  (Continued).   Listing  of  the  IIP9830  Index  of 
Refraction  (n-1)  Program. 


513 
520 
539 
540 
550 
560 
570 
530 
590 

600 

610 
620 
6  3  0 
640 
650 
6  6  0 
670 
6  3  0 
690 
691 
692 
6  9  3 

700 

705 
710 
720 
730 
740 
770 
790 
3  0  0 
310 
320 
350 
3  3  0 
9  0  0 
910 
920 
930 
960 
930 
9  9  0 

100 

102 
103 
104 
105 
106 
109 


i'  - 
GC 
GC 
GC 
GC 
GC 
EC 
EC 
EC 
EC 
EC 
EC 
EC 
EC 
EC 
EC 
EC 
EC 

r-.  i 

CI 
H  = 


0 

1  3=1 

2  3=GC4  3  =  GC7  3  =  GC  12  3  =  1 

3  3=GC6  3  =  GC  10  3  =  36 
5  3=GC3  3='GC  11  3=73 
9]=34 

1   "J  =  0 

2  3=134674,,  61 
=150197.31 

=153950. 1 
=161374.4 
=163130. 29 

=  16650  0 
=166562 
3=167312.5 
3=16  8 0  0  0 


3  3 

4  3 
53 

5  3 
73 
3  3 
9  3 

10 

11 
12 


4. 


3=163750 
.6171E-05 
. 997925E+10 
1356E-15 


C2=C1*H 

C=lE-27 

Y  =  Y+1 

FOR  ,1=1  TO  41 

READ    TCJ3jNCJj  1  3jNCJj2 

NEXT    J 

FOR    J=l    TO    41 

NEXT    J 

FOR    J=l    TO    41 

Q=NC  Jj2  3*GC  1  ]/<NC  Jj  1  3* 

Q=LOGQ 

KC  J3  =  -((EC2  3-EC  1  3>*C2) 

NEXT    J 

FOR    J=l    TO    41 

FOR    K=4    TO    12 

Ql  =  ((ECK3-EC3  3>*C2).--a< 

Q1=EXP(-Q1> 

RCJjK3=GCK3*Q1/GC33 

NEXT  K 

NEXT  J 

FOR  ,1=1  TO  41 
0  3=1 

0    FOR    K=4    TO    12 
0    S=S+ACJjK3 
O    NEXT    K 
O    SCJ3  =  3 

0    NC  Jj3  3  =  ZC  J3/SC  J3 
O    NEXT    J 


3jZCJ3 


:c2 


3) 

:c  J3> 
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TABL"  IX  (Continued).   Listing  of  the  HP9830  Index  of 
Refraction  (n-1)  Program. 


1  100 

FOR  ,1=1 

TO  41 

1  120 

FOR  K  =  4 

TO  12 

1 130 

HCJjK]=1 

«J*33*fl[J»K 

1 160 

NEXT  K 

1 170 

NEXT  J 

1210 

MC  1  ]  =  55i 

~\ 

1220 

MC2]  =  -3" 

7444 

1230 

MC3]  =  -5' 

3674 

1240 

MC4  3=-5i 

-1 7  '"' '"' 

1250 

MC  5  3=31 

324 

1260 

NC6  1  =  -6 

1753 

1270 

MC7]  =  -5 

4913 

1230 

MC3  3  =  27 

317 

1290 

N[9]  =  0 

1300 

MC  10]  =  - 

32435 

1 3 1 0 

MC  1 1 ]  =  4 

3799 

1 3  2  0 

MC  12  1  =  3 

4233 

1330 

FOR  J=l 

TO  41 

1340 

R  =  0 

1350 

FOR  K=l 

TO  12 

1 3  6  0 

RC  JjK3= 

MCK3*N[J»K3* 

1370 

R  =  R  +  R[,I 

t  K  ] 

1400 

NEXT  K 

1412 

IF  Y=l 

THEN  1413 

1414 

HCJ]  =  R 

1416 

GOTO  14 

20 

1-413 

LCJ]  =  R 

1420 

PRINT  " 

N  -  1  IS"R"H 

1440 

NEXT  J 

1450 

PRINT 

1460 

PRINT  " 

HIGH 

1465 

PRINT 

1463 

IF  Y  >  = 

2  THEN  1500 

1470 

GOTO  70 

c 
,1 

1  5  0  0 

SCALE  O 

p 2000? -2E-05 

1510 

X  FIX  IS  - 

2E-05? 100 

1520 

Y fix  is  o 

h 1E-05 

1530 

FOR  ,1  =  1 

TO  41 

1540 

PLOT  TC 

J3»LC  J] 

1550 

NEXT  J 

1560 

PEN 

1570 

FOR  ,1=1 

TO  41 

1530 

PLOT  TC 

J3jHCJ3 

1590 

NEXT  J 

1  6  0  0 

PEN 

1610 

STOP 

TINE 


:  J  3 "MSEC 

ZR    llfl/CI'1t2 


6E-05 
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TA3L2  IX  (Continued).   Listing  of  the  HP9830  Index  o 
Refraction  (n-1)  Progran. 


c 


5090  PRINT  " 

5010  DATA  0) 

5029  DATA  50 

5030  DATA  10 
5040  DATA  15 
5050  DATA  20 
5060  DATA  25 
5070  DATA  36 

c,  d  O  i%  fiuTQ  "'  =; 


50? 


5 1 89  DATA  850) 9 . 7 1 SE+ 1 9  ? 1 . 1 73E  + 1 4?  6. 755E  +  1  3 

5190  DATA  900?  9. 71SE+19? 1 „ 174E+14) 6, 756E+13 
5200  DATA  950j9.716E+19j I . 174E+14) 6. 757E+13 
5210  DATA  1000? 9.71SE+19) 1 . 175E+14? 6. 757E+1 
5220  DATA  1950? 9. 71SE+19) 1  - 175E+14?  6. 758E+1 
5230  DATA  1 1Q0> 9. 716E+19 j 1 . 175E+14? 6. 758E+1 
5240  DATA  1 15© j 9. 716E+19) 1 . 17SE+14) 6. 759E+1 
5250  DATA  120O?  9.716E+19?  1  .  176E+14)  6.  76E+1S 
5260  DATA  1 250? 9. 7 1 6E+ 1 9 ? 1 . 177E+ 1 4? 6 . 76E- 
5270  DATA  1300? 9. 7 16E+ 19? 1 . 173E+ 14? 6. 67  1 E 
5230  DATA  1350? 9 . 716E+1 9 > 1 . 17SE+1 4 ? 6 . 7621 
5290  DATA  1400 ? 9. 716E+ 19? 1  - 179E+14? 6 . 763E 
5300  DATA  1450?9.716E+19? 1 . 13E+14? 6. 765E- 
5310  DATA  1500? 9. 716E+19? 1 . 181E+14? 6. 766! 
5320  DATA  1550? 9 . 716E-' 
5330  DATA  1600?9.716E 

=;  ■"!  A  u  n  Q  T  H  IkSu.Q  7  1  hF 


5340  DATA  165i 
5350  DATA  17Q| 
5360  DATA 
5370  DATA 
5330  DATA 
5390  DATA 
5400  DATA 
5410  DATA 
5420  END 
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TABL2   X.      Output   of   the   HP9330   for  n-1    at  Low  Power 


LOU    POWER    lOfl/CII' 


H  -  1 

IS 

5. 42153E-05 

AT 

TIME 

8 

NS 

;EC 

N  -  1 

IS 

5.42126E-05 

A  T 

TIME 

50 

N:; 

:EC 

H  -  1 

IS 

5. 42097E-05 

RT 

TIME 

100 

m 

;EC 

H  -  1 

IS 

5, 42878E-85 

AT 

TINE 

150 

H:i 

>EC 

H  -  1 

T  '«' 

5.42867E-05 

AT 

TIME 

200 

NS 

SEC 

M  -  1 

IS 

5.428  6 U  E - 0 5 

AT 

TIME 

2  5  0 

l-l:- 

5  EC 

H  -  1 

IS 

5. 42856E-05 

AT 

TIME 

2  U  @ 

H:l 

5EC 

N  -  1 

T  C* 

5. 42053E-05 

RT 

TIME 

350 

HS 

;;ec 

H  -  1 

T  C* 

5.42052E-05 

RT 

TIME 

4  0  0 

N' 

5  EC 

N  -  1 

IS 

5. 42051E-85 

RT 

TIME 

459 

N: 

::EC 

N  -  1 

IS 

5. 42050E-05 

RT 

TIME 

5  0  0 

N: 

LEG 

N  -  1 

IS 

5. 42050E-85 

RT 

TIME 

C"  l~  |"i 

N: 

3  EC 

N  -  1 

T  '«* 

5.42050E- 0  5 

Hi" 

TIME 

6  8  0 

!■  I: 

:.  EC 

N  -  1 

IS 

5.42058E-05 

AT 

TIME 

650 

N: 

3  EC 

N  -  1 

IS 

5. 42058E-05 

RT 

TIME 

7  0  0 

H: 

'i  EC 

N  -  1 

T  ■— ' 

5. 42049E-85 

AT 

TIME 

-?  cr  |"i 

N: 

SEC 

N  -  1 

IS 

5. 42049E-05 

RT 

TIME 

:  ■;:  0  0 

H: 

SEC 

N  -  1 

T  C ■ 

5. 42049E-85 

AT 

TIME 

350 

M: 

SEC 

H  -  1 

I  s 

5. 42049E-05 

AT 

TIME 

1 0  8 

N: 

SEC 

N  -  1 

IS 

5.42049E-85 

AT 

TIME 

9  5  0 

H 

SEC 

N  -  1 

T  '-' 

5 .  4  2  8  4  9  E  -  0  5 

R  T 

TIME 

1  0  0  0 

\\'. 

;EC 

H  -  1 

IS 

5  .  4  2  0  4  9  E  -  0  5 

RT 

TIME 

1050 

NS 

SEC 

N  -  1 

T  '— ' 

5 . 4204  9 E - 0  5 

RT 

TIME 

1  100 

MS 

SEC 

N  -  1 

T  '"* 

5 .  4  2  0  4  9  E  -  0  5 

RT 

TIME 

1150 

NS 

SEC 

H  -  i 

IS 

5.42049E-05 

RT 

TIME 

1288 

N'c 

SEC 

H  -  1 

IS 

5.42049E-05 

RT 

TIME 

1250 

NS 

SEC 

N  -  1 

T  '— ' 

5. 42049E-05 

RT 

TIME 

1  3  0  0 

m 

3  EC 

H  -  1 

IS 

5.42049E-05 

AT 

TIME 

1350 

M' 

SEi. 

N  -  1 

T  '— ' 

5. 42049E-05 

AT 

TIME 

1400 

N: 

SEC 

N  -  1 

T  C* 

5. 42049E-05 

AT 

TIME 

1450 

N: 

SEC 

N  -  1 

IS 

5. 42049E-05 

AT 

TIME 

1500 

N: 

SEC 

N  -  1 

IS 

5.42849E-05 

AT 

TIME 

1550 

N: 

SEC 

H  -  1 

IS 

5. 42049E-U5 

AT 

TIME 

1608 

N: 

SEi 

N  -  1 

T  '— " 

5. 42049E-05 

RT 

TIME 

1650 

H: 

SEl 

N  -  1 

T  '— * 

5. 42849E-85 

AT 

TIME 

1700 

H: 

SEC 

H  -  1 

IS 

5.4204  9  E - 0  5 

RT 

TIME 

1750 

N: 

SEL 

N  -  1 

IS 

5.42849E-05 

RT 

TIME 

1  8  0  0 

H: 

SEi 

N  -  1 

IS 

5.42049E-05 

AT 

TIME 

1350 

H: 

SEi' 

N  -  i 

T  '— ' 

5. 42O43E-05 

AT 

TIME 

1  9  0  0 

N: 

SEi 

N  -  1 

IS 

5. 42048E-85 

AT 

TIME 

1958 

N: 

SEi 

N  -  1 

IS 

5. 42048E-05 

AT 

TIME 

2  8  0  0 

H 

3  Ei 
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TABLZ  XI.   Output  of  the  HP9330  for  n-1  at  High  Power 


HIGH  POWER  llfl. 


t2 


H  -  1 

IS 

5.42153E-05 

RT 

TIME 

0 

HS 

>E 

j 

H  -  1 

IS 

5.42123E-05 

RT 

TIME 

50 

Hi 

*E 

;: 

N  -  1 

IS 

5.42091E-05 

ht 

TIME 

100 

N:! 

5E 

"• 

N  -  1 

IS 

5.42071E-05 

RT 

TIME 

1 5  0 

N:- 

5E 

H  -  1 

I  s 

5 .  4  2  0  5  8  E  -  0  5 

H  T 

TIME 

200 

H:; 

5E 

H  -  1 

T  C" 

5.42051E-05 

RT 

TIME 

250 

H: 

;;E 

z 

N  -  1 

IS 

5. 42046E-05 

HT 

TIME 

3  0  0 

H: 

5E 

2 

N  -  1 

IS 

5. 42043E-U5 

RT 

TIME 

350 

H: 

5E 

2 

H  -  1 

T  '-•' 

5. 42042E-05 

RT 

1  1ME 

400 

M: 

:>t 

2 

H  ~  1 

T  '-' 

5.42041E-05 

RT 

TIME 

450 

N: 

3E 

2 

N  -  1 

I  s 

5  .  4  2  0  4  0  E  -  0  5 

RT 

TIME 

5  0  0 

N: 

3E 

2 

N  -  1 

IS 

5. 42040E-05 

RT 

TIME 

550 

N: 

E;E 

2 

H  -  1 

I  s 

5. 42039E-05 

RT 

TIME 

6  0  0 

M: 

:>  L 

2 

H  -  1 

T  C' 

5.42039E-05 

RT 

TIME 

650 

N: 

3E 

2 

H  -  1 

I  s 

5 .  4  2  9  3  9  E  -  0  5 

RT 

TIME 

7  0  0 

H: 

3E 

c 

N  -  1 

T  C* 

5 .  4  2  0  3  9  E  -  9  5 

R  1" 

TIME 

-?  cr  i-j 

i  -.jU 

N: 

5E 

c 

N  -  1 

IS 

5. 42039E-05 

RT 

TIME 

8  0  0 

N: 

3E 

c 

H  -  1 

IS 

5.420  3 9 E - 0  5 

RT 

TIME 

i~l  C  i~j 

O  J  <0 

N 

5E 

c 

N  -  1 

T  '"' 

5. 42039E-05 

RT 

TIME 

9  0  0 

N 

BE 

c 

N  -  1 

IS 

5 . 4203 9 E - 0  5 

R  T 

TIME 

'■i  5  0 

H 

3E 

c 

H  -  1 

IS 

5. 42939E-05 

RT 

1  1ME 

1  0  @  0 

NShU 

N  -  1 

IS 

5 .  4  2  0  3  9  E  -  0  5 

HT 

TIME 

1  0  5  0 

MSEC 

H  -  1 

IS 

5. 42039E-05 

RT 

TIME 

1  1  0  0 

MSEC 

N  -  1 

T  '— ' 

5. 42039E-05 

RT 

TIME 

1150 

MSEC 

N  •-  1 

T  C* 

5. 42039E-05 

RT 

TIME 

1  2  0  0 

MSEC 

N  -  1 

IS 

5 .  4  2  0  3  9  E  -  0  5 

RT 

TIME 

1250 

MSEC 

H  -  i 

IS 

5 .  4  2  0  3  9  E  -  0  5 

RT 

TIME 

1  3  i  j  0 

MSEC 

N  -  1 

T  C* 

5. 42039E-05 

RT 

TIME 

1350 

MSEC 

N  -  1 

T  '-' 

5. 42039E-05 

RT 

TIME 

1  4  0  0 

MSEC 

H  -  1 

IS 

5. 42033E-U5 

RT 

TIME 

1450 

MSEC 

N  -  1 

T  '-* 

5. 42033E-05 

RT 

TIME 

1  5  0  0 

MSEC 

N  -  1 

T  C- 

5.42038E-05 

RT 

TIME 

1550 

MSEC 

N  -  1 

T  '«' 

5. 42033E-05 

RT 

TIME 

1  6  0  0 

MSEC 

H  -  1 

T  '-' 

5. 42033E-05 

RT 

TIME 

1650 

MSEC 

N  -  1 

IS 

5.42033E-05 

RT 

TIME 

1700 

MSEC 

H  -  1 

IS 

5. 42038E-05 

RT 

TIME 

1750 

MSEC 

N  -  1 

T  '- 

5.4203SE-05 

RT 

TIME 

1 8  0  0 

MSEC 

H  -  i 

T  C' 

5.42033E-05 

RT 

TIME 

1850 

MSEC 

N  -  1 

T  C' 

5.42037E-05 

RT 

TIME 

1  ^  0  0 

MSEC 

N  -  1 

T  '- 

5.42037E-05 

A  T 

TIME 

1950 

MSEC 

H  -  1 

I S 

5.42037E-05 

RT 

TIME 

2  0  0  0 

MSEC 
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VII.   CALCULATION  OF  PHASE  SHIFT 

Beam  quality  as  previously  mentioned  depends  on  favor- 
able illumination  of  the  exit  aperture;  favorable  illumina- 
tion implies  near  constant  phase.   Therefore,  calculation 
of  the  phase  shift  is  important.   To  calculate  the  phase 
shift  consider  a  laser  cavity  of  length  L  as  shown  in  Figure 
10. 


CAVITY 


LASER  OUTPUT 


Figure  10.   Geometry  of  Cavity 

Other  symbols  essential  for  the  calculation  are 

N  -  number  of  waves       A  -  wavelength 

1  -  denotes  low  power     n  -  index  of  refraction 

2  -  denotes  high  power    v  -  frequency  of  the  laser 

The  number  of  waves  in  length  L  is  given  by 


N.  -  L/X. 


(59) 
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In  equation  (59)  i  =  1,2.   Wavelength  is  related  to  index 
of  refraction  by 


X.    =  v./v  =  -£-  (60) 

i    l     n .  v  v   ' 


Consider  two  regions  in  the  cavity  which  are  pumped  at 
different  power  inputs.   The  difference  in  the  number  of 
waves  in  the  two  regions  is 

N2  "  ?'l  "  L<4  "  *?  (61) 

Substituting  equation  (60)  into  equation  (61)  yields 

n0v   n, v    T 
U2  -  ^  =  L(- L.)  =  If    („2-ni)  (62) 

Multiplying  the  right  hand  side  of  equation  (62)  above  and 
below  by  n-.  produces 

Lv  n,   n9-n, 

N?  "  Ni  =  —T1    C-£— -)  ^3) 

I  1     c     n^ 

Multiplying  equation  (63)  by  1/L  and  substituting  equation 
(50)  yields 


N2  -  N1    ,   n«-n, 

-i-r-^  -  ^  Hq-1'  <64> 


The  term  on  the  left  hand  side  is  the  phase  shift  per  unit 
length  of  cavity.  The  phase  shift  at  1000  nanoseconds  for 
a  nonresonance  wavelength  is  now  calculated. 

X   =  353  nm  =  353.0E-7  cm  (Figure  8) 
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From  Tables  X  and  XI  the  values  of  n  are 


n£  =  1  +  5.42039S-05 


n,  -  1  +  5.42049S-05 


Substituting  into  equation  (64) 


N2"N1       1     ,1.00005420 39 -1.0000542049^     9  s„7„  _    -1 
L   '  353.0E-7  k     1.0000542049         ;    -z.  ojz/£,-:>  cm 


The  phase  shift  at  1000  nanoseconds  for  resonance  is  calcu- 
lated below. 

X   =    351.6E-7 

The  values  of  n  are 


n2  =  1  +  (-3.71283E-06) 


nl   =  1  +  (1.6602E-06) 
Substituting  into  equation  (64) 


N2~^l       1     ,.999996287170-1. 000005602  v     «  Q«-QT?  ,    -1 
L   "  351.6E-7  K  1.000006602         ;    -a.bO^-J  cm 


The  Strehl  ratio,  which  is  discussed  by  Born  and  Wolfe 

[21] ,  provides  a  criterion  for  beam  quality  in  terms  of 
phase  distortion.   The  ratio  is 

I  -  I   exp  [-(2tt6/A)2]  (65) 
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where  I  is  beam  intensity  in  the  far  field  due  to  a  phase 
distortion  of  6/Xo  Based  on  equation  (65),  the  criterion 
of 

6/X  <  0.1  (66) 

is  frequently  invoked  for  satisfactory  beam  quality.   Com- 
paring equations  (65)  and  (61)  one  notes  that 

N2  -  tf-  -  6/X  (67) 

The  maximum  length  of  the  laser  cavity,  L,  for  good 
beam  quality  is 

L  =  ^_   =   _J — !_  (68) 

n2-n1      Un-U 
{-V^-}     (21) 

Substitution  of  numerical  values  yields  the  following 
results : 

-  maximum  L  for  a  laser  operating  at  a  nonresonant  wave' 
length  is  35  meters . 

-  maximum  L  for  a  laser  operating  at  a  resonant  wave- 
length is  0.1116  meters. 

For  a  nonresonant  laser  wavelength,  a  laser  cavity  with 
length,  L,  less  than  35  meters  will  have  6/X    less  than  0.1. 
However,  if  the  laser  wavelength  is  in  resonance  with  neon, 
a  length  of  laser  cavity  of  only  0.1116  meters  will  cause 
5/ A  =  0.1.   Recall  the  power  into  the  laser  cavity  was 
varied  by  10  percent. 


66 


VIII.   RESULTS 

For  a  power  variation  of  10  percent  the  phase  shift  at 
nonresonance  is  small  (N2  -  N-,  =  -0.0028  for  a  100  cm 
cavity)  and  is  dominated  by  the  neon  ground  state  whose 
population  is  approximately  one  million  times  larger  than 
the  other  species.   However,  in  the  resonance  case,  if  the 
assumed  values  of  R*  are  correct,  the  index  of  refraction 
can  be  greatly  magnified,  resulting  in  a  phase  shift  of 
-8.959S-3  cm   as  compared  to  -2.3327E-5  cm   for  nonreso- 
nance . 

It  should  be  noted  that  the  populations  of  the  neon 
excited  states  were  calculated  using  an  arbitrary  choice  to 
determine  the  rate  constants.   However,  the  population  of 
Ne  is  approximately  six  orders  of  magnitude  larger  than 
the  neon  excited  states  and  thus  dominates.   As  long  as  the 
laser  is  operated  away  from  the  3s  -*■  4p  resonant  wavelengths, 
the  phase  shift  will  be  negligibly  small  resulting  in 
satisfactory  beam  quality. 
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APPENDIX  B 
CALCTJLATION  OF  RATE  CONSTANT  k5 

Consider  a  thin  layer  of  active  optical  medium,  with 

gain  a,  cm   ,  in  the  laser  cavity  of  length  I,    cm.   The 

2 
change  in  intensity  I,  watts/ cm  ,  within  the  thin  layer  is 


The  number  of  photons  per  unit  time  entering  the  layer  of 
active  optical  medium  is 


I.   A 

N.   =  tn  (79) 

in    nv  v  y 


A  corresponding  term  for  the  number  of  photons  leaving  per 
unit  time  is 


I   .  A 

N    =   out  (71) 

out     hv 


The  rate  of  production  of  photons  within  the  thin  layer  is 

•  • 

the    difference   between   N  and  N.      which   is 

out      in 

( rate  of  pro-  )  •      A 

duction      =  N   ,.  -  N.   =  A  (I   _  -  I.  )        (72) 
}  -  ,  .      (     out     in   hv  v  out    m'        v 
( of  photons   } 

Equation  (72)  neglects  spontaneous  emission.   Combining 
equations  (69)  and  (72)  yields 
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production  rate  ■  r~  (73) 


From  equation  (19)  the  rate  of  production  of  photons  by- 
stimulated  emission  is 

production  rate  -  k^thv]  [KrF*]A£,  (74) 

where  the  symbol  [hv]  means  number  of  photons  per  unit  volume 
The  gain  coefficient  also  is  equal  to 

a  =  a[KrF*]  (75) 

where  a  is  the  optical  cross  section  for  stimulated  emission 
at  the  frequency  v.   The  number  density  of  photons  [hv]  is 
related  to  intensity  by 

[hv]  -  -g_  (76) 

where  c  is  the  velocity  of  the  photons;  c,  obviously,  is  the 
speed  of  an  electromagnetic  wave. 

Combining  equations  (73)  to  (76)  gives 

Ic5  =  ca  (77) 

3 

The  dimensions  of  k-  are  cm  /sec;  the  dimension  of  a  is  cm, 

2 
and  c  is  cm  /sec.   From  reference  [2]  the  value  of  a   is 

2 
5. 011-16  cm  ,  and  c  is  3.0E10  cm/sec.   Substituting  these 

values  into  equation  (77)  yields 
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k5  =  1.5E-05  cm3/sec 


This  is  the  value  for  k_  used  in  the  computer  program 
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APPENDIX  C 

CALCULATION  OF  PHOTON  CONCENTRATION 

The  concentration  of  photons  is  given  by  equation  (76) 
as 

[hv]  =  I/chv 


2 
Assuming,  I   =    10   kcu/cm      (saturation  value) 

c   =    3.0E10    cm/sec 

h   =    6.6256E-34  J-sec 

XQ   =    243.0E-9   m 


the  concentration  of  photons  can  be  calculated. 


r,  -,   .  10.QE3 

LhvJ  -  *  riF3"" 

(3.0E10)  (6.6256E-34)  (2te.VS.P 


[hv]    =   4.159Ell/crn3    . 
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